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Abstract 

A new approach to the theory of polynomial solutions of g-difference equations is 
proposed. The approach is based on the representation theory of simple Lie algebras 
Q and their q - deformations and is presented here for U q {sl(n)). First a g-difference 
realization of U q {sl(n)) in terms of n(n — l)/2 commuting variables and depending on 
n — 1 complex representation parameters r^, is constructed. From this realization lowest 
weight modules (LWM) are obtained which are studied in detail for the case n = 3 (the 
well known n = 2 case is also recovered). All reducible LWM are found and the polynomial 
bases of their invariant irreducible subrepresentations are explicitly given. This also gives 
a classification of the quasi-exactly solvable operators in the present setting. The invariant 
subspaces are obtained as solutions of certain invariant g-difference equations, i.e., these are 
kernels of invariant g-difference operators, which are also explicitly given. Such operators 
were not used until now in the theory of polynomial solutions. Finally the states in all 
subrepresentations are depicted graphically via the so called Newton diagrams. 
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1. Introduction. 

There has been impressive progress in the last few years toward solving a fundamental 
problem posed by Bochner [1]: characterize differential operators possessing orthogonal 
polynomial eigenfunctions. This progress has taken many guises corresponding to several 
very different areas of research. In physics, the motivation has been to obtain solvable 
interaction Hamiltonians for bound state [2] and scattering problems [3] of the Schrodinger 
equation by using Lie-algebraic techniques (defined precisely below). Such techniques 
have re-appeared recently in the discovery by Turbiner, [4], Shifman [5], and others, of a 
new class of physically significant spectral problems, the so-called "quasi-exactly solvable 
problems", in which part of the (bound state) spectrum can be found algebraically. 
The basic structure of these related approaches has been recently clarified by the more 
mathematically oriented work of Kamran and Olver [6], [7], [8]. In classical applications of 
symmetry in quantum physics, using Lie group-theoretic methods, the Lie group appears as 
a symmetry group leaving the Hamiltonian invariant. A more general point-of-view arose 
in the late sixties with the introduction of spectrum generating groups [9] (and algebras) 
in which the symmetry group no longer necessarily left the Hamiltonian invariant but 
still retained some symmetry structure in the sense that the spectrum was "split but 
not mixed" in Gell-Mann's phrase [10]. In this more general approach, a second-order 
differential operator - - the Hamiltonian - - admits a spectrum generating symmetry if 
the Hamiltonian can be written as a quadratic in the generators (first-order differential 
operators) of a finite dimensional Lie algebra. The surprisingly successful interacting boson 
model of Iachello and Arima [2] is of this form with a compact Lie group symmetry (SU(6)); 
non-compact Lie group symmetries (SX(3, -R), Sp(6, M), . . .) have also enjoyed favor [11]. 

Following Kamran and Olver [7] we define a differential operator to be Lie-algebraic if 
it can be written as a quadratic combination of first-order differential operators (vector 
fields) which generate a finite-dimensional Lie algebra. A finite-dimensional Lie algebra is 
said to be [7] quasi-exactly solvable if it has a finite dimensional representation on some 
subspace of the space of smooth functions. In this terminology, a quasi-exactly solvable 
differential operator is a Lie-algebraic differential operator corresponding to a quasi-exactly 
solvable Lie algebra (the "hidden symmetry algebra"). As discussed by Shifman and 
Turbiner [5], the spectrum associated with the finite-dimensional representation spaces 
can be computed algebraically. 

As one might expect, the one-dimensional case has been solved completely. It is useful 
to review these results here, since this case shows most clearly the essential underlying 
ideas: Theorem [7] (Gonzalez-Lopez, Kamran, Olver): Every finite- dimensional quasi- 
exactly solvable Lie algebra of first order differential operators in one (real or complex) 







variable is locally equivalent to a subalgebra of one of the Lie algebras 

( d d 2 d 1 
g r = Span z — , z — - rz, 1 > , (1) 

[ dz dz dz J 

where r £ is a non-negative integer. The associated representation module consists 

of polynomials of degree r. 

A cohomological interpretation of the integer r has been given by Gonzalez-Lopez et 
al [8]. (The fact that r must be integral has been called "quantized cohomology", and 
these authors [8] suggest the underlying reason as a problem for further investigation.) 

The words "locally equivalent" in this theorem has a precise and important meaning: 
Two differential operators are said to be locally equivalent if they can be mapped into 
each other by a combination of change of independent variable: x = f(z), and a (not 
necessarily unitary) "gauge transformation": O — > e 9 ^ Oe~ 9 ^ z \ where O is a differential 
operator. 

As an abstract Lie algebra g r in the result above is a central extension of the subal- 
gebra: 

f d r 2 d <i 1 

K " \ Jo = Z dz--2 ' J+=TZ - Z Tz> J - = Tz>) ■ (2) 
It is sufficient, as Kamran and Olver remark, to concentrate attention - - without loss of 
generality - - on the algebra h r rather than on the full algebra g r , since any Lie algebraic 
operator of g r is automatically a Lie-algebraic operator of h r . 

The Lie algebra given by (2) is very familiar to physicists: it is just the Lie algebra 
with the commutation relations: 

[Jo, J±] = ±J± , [J+,J-] = 2J , (3) 

which generate the quantal angular momentum group, SU(2). The realization of this 
algebra (which can be expressed in terms of one pair (a + ,a) of boson operators), is 
probably less generally familiar but equally important: 

Remark: The realization, (2), is the simplest example of a general algebraic technique 
for constructing all_ unitary irreps of all compact simple Lie groups, the algebraic Borel- 
Weil construction. 

The importance of this remark belies its ease of proof [12], [13], [14]. It is easily seen, 
for example, that the "quantization of cohomology" found in [8] is a direct consequence 
of the Borel-Weil construction (compactness) and as such appears to require no separate 
proof. 

A further consequence of our remark is that the algebraic Borel-Weil construction im- 
plies a generalization of Lie algebraic quasi-exactly solvable Hamiltonians to encompass all 
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simple Lie algebras. Whether or not this generalization encompasses all such structures, 
in the sense of the GKO theorem, is a question for further investigation. 

The present paper is concerned with an application of our remark to the simple Lie 
algebra Ag, though it is clear, that the ideas are easily extendable for any simple Lie 
algebras Q. As a rank £ Lie algebra, the Cartan toroid must involve £ integers (rather 
than one as in the one-dimensional {A\) algebra discussed above). Correspondingly, the 
algebraic Borel-Weil construction for At must involve £(£ + l)/2 complex variables 
(commuting bosons, which parametrize a Borel subgroup of SL(£ + 1) [15]), thus going, 
even for £ = 2, beyond the recent investigations of two variables [8]. 

We have chosen to discuss our extension of quasi-exactly solvable systems to more 
than two variables by using the more general quantum group U q {sl(n)) instead of s/(n), 
where q is a generic parameter. This is technically a much more difficult construction, 
but it is at the same time much more general. Results for the Lie algebra sl(n) are easily 
obtained by specializing our results to q = 1. 

Note also that all our considerations are valid also for U q (sl(n } 1R)) if we restrict to 
real variables, real representation parameters and real q. Nothing is lost from the structure 
of the representations since all reducibilities occur when the representation parameters or 
their combinations are nonnegative integers. 

We now give a summary of our approach and the results of the paper. The approach 
has several ingredients which are in fact available via the representation theory of simple 
Lie algebras Q and their q - deformations. 

The first ingredient is the construction of a q - difference realization of the repres- 
entations of U q {Q). For U q {sl(n)) we use the construction of [15] for U q {su(n)). The 
additional input w.r.t. [15] is to assign explicitly the eigenvalues of the Cartan gener- 
ators. These are the representation parameters of our representations. Moreover, these 
eigenvalues are taken to be arbitrary complex numbers in order to obtain the most gen- 
eral representations. We also express all quantitites in terms of n(n — l)/2 commuting 
variables mentioned above (n = £ + 1), and q - difference and number operators in these 
variables. 

Applying this realization to the constant function (independent of any of the variables) 
we obtain in fact realizations of lowest weight modules of U q {sl(n)). In the present paper 
we consider the case when q is not a nontrivial root of 1. Matters are arranged so that for 
generic representation parameters r^, z = 1, . . . n — 1, these representations are irredu- 
cible and infinite-dimensional. If some representation parameters (or certain combinations 
thereof) are nonnegative integers the representation is reducible. The invariant irreducible 
subspace, or subrepresentation, is finite-dimensional iff all representation parameters are 







nonnegative integers: V{ G -Zv+. In all cases the subrepresentations have polynomial bases. 

Already at this point we can make a classification of all quasi-exactly solvable op- 
erators, which is equivalent to classifying the subrepresentations. This classification will 
depend entirely on the values of the representation parameters. Such a classification is 
also equivalent to the multiplet classification approach developed in [16]. In particular, 
for specific values of the representation parameters r^, namely, only when some of them 
are zero, the invariant subspaces will depend on less variables; however, they can still be 
infinite-dimensional. As in the case of functions in one and two variables [17] all this can 
be visualized using the relevant Newton diagrams [18] which for the finite-dimensional 
representations are certain convex polyhedra in IE!^ n 1 ^ 2 . Moreover, it is natural from 
our point of view to introduce also infinite Newton diagrams corresponding to infinite- 
dimensional nontrivial subrepresentations. 

We also discuss the g-difference operators which determine the g-difference equations 
whose solutions are the invariant subspaces. These are operators which arise when our 
representations are reducible and which intertwine these representations. We should stress 
that such operators, though widely used in many branches of mathematics and physics, are 
not used in the literature mentioned above. In general, such operators correspond to each 
positive root of the root system of Q arising whenever the corresponding representation 
parameters are nonnegative integers. This situation is well understood in the q = 1 case 
[19], and the algebraic part of it is developed also for the U q {Q) case [20]. 

The above approach is implemented fully in the present paper to the case U q (sl(3)) (cf. 
Sections 3.3-3.6, 4.), while the general case is postponed to a sequel of the paper. (In 
passing, we also recover the well known U q {s!(2)) case (cf. Section 3.2.)) For U q (sl(3)) we 
give the basis of all invariant subspaces. The explicit polynomial basis (in three variables) 
needs hypergeometric functions for compact notation. Here there are three possible in- 
variant operators of which at most two exist or are relevant for a particular invariant 
subspace. There are two generic situations, one of which contains, as a speaial case, the 
finite dimensional representations. The Newton diagrams are given in Section 4. 
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2. Procedure for the construction of the representations. 

The procedure is iterative. In fact, we have to use also U q {gl{n)). Let us introduce 
first some notation. The basic g-number notation is: 

a/2 _ -a/2 

and it will be used also for diagonal operators H replacing a in (4). Following [15] our 
representations will be given in terms of n(n — l)/2 variables. For our purposes we denote 
these variables by zf , 2<k<n,l<i<k — 1. Next we introduce the number operator 
for the coordinate z\ , i.e., zj 1 = Sm^SijzJ 1 and the q - difference operators D\ , 
which admit a general definition on a larger domain than polynomials, but on polynomials 
are well defined as follows: 

D * = ^[iVf] • (5) 
z i 

Further we note that the representations of U q {sl(n)) will be characterized by n — 1 com- 
plex parameters r& G (F, 1 < k < n — 1. 

We rewrite formulae (5.3), (6.10) and (6.22) from [15] in the following way : 

T n (E tJ ) = Tn^Etj) q^- N P + qkTn-^Ejj-Eu) ^ ? i<j <n , (6a) 
T n (E ZJ ) = T n -i(Eij) qW-ND + qFn-^-Ess) z? jjn ? n>l> ^ { Q b) 
T n (Eu) = Tn-^Eu) + iVf , z < n , (6c) 

r n (^ nn ) = r n _!(E nn ) - , ( 6c/ ) 

i=l 

r n (^ nJ ) = 9 ^(E]:i r »-i( s »)-E J ": i 1 + i r »-i( s »)) ^ 5 2<n? ( 6e ) 



r n (^ m ) =g a » < [r^^)-^-!^)- Viv, 



/ y - A: 
fc=l 



n-1 

i=i 
i 5^1 



(6/) 



where 
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( + )( l(N n + N n )) (^H (E ''<*<■> if ''<■> ^fc) 

+ ( + for z < j , 

O") = 1 ~ for z > j , 
I for i = j . 



Y n -i(Eij) are defined at the previous step, except Y n -\(E nn ) which adds the representa- 
tion parameter r n -\ and is given by: 

n-l 

Tn-^Enn) = = + T . (8) 

A;=0 

The parameter ro represents the center of U q {gl{n)) and is decoupled later. Note that 
(6a — <f), (6e), (6/), respectively, correspond to (5.3a-d), (6.10), (6.22)) of [15], respect- 
ively. The additional input with respect to [15] is : 1) notation - we put indices on Y 
corresponding to the case we consider - thus this is made an iterative procedure; 2) we 
give values to the Cartan generators which values are consistent with previous knowledge 
from representation theory; 3) we introduce q - difference operators D™ to replace z™\ 
4) we have done also something artificial by using Y n -\{E nn ) - this is given by the 'total 
number' r n_1 (which for finite dimensional representations is the number of boxes of the 
Young tableaux minus n — 1) plus the number ro representing the center of U q {gl{n)). 

Denoting Z n = ^T=i we have: 

n 

Y n (Z n ) = ]T T n (Eu) = r n ~ x + r + Y n - X {Z n - x ) = 

(9) 

n-l n-l v ' 

= r l + nro = (n — i)ri . 

i=l i=0 

Thus, as expected, Y n (Z n ) is central. Then the generators H™ = Y n (Ea — Ei+i^+i), 
1 < i < n, Y n (Eij) } i ^ j, form a q - difference operator realization of U q {sl(n)). 

It is straightforward to obtain the explicit expressions for Y n (Eij). In particular, we 

have 

n — i — 1 i — 1 i — 1 

rn(£«) = E N *~ J - T. N i + E r > ' ( 10 ) 

j = j=l j = 

with the usual convention that a sum is zero if the upper limit is smaller than the lower 
limit. From this we obtain the expressions for the Cartan generators H™ (as defined above): 

n—i—2 i—1 

H? = 2Nl +1 - n + ( N r 3 -N^) + E(iVj +1 -iV;) z<n. (11) 

j = j=l 

Let us illustrate things for n = 2, 3. 

For n = 2 we have (we use only (6e, /), (11)) : 

r 2 (£ 12 ) = zl[ ri -Nl] = x[r-N x ], (12a) 
Y 2 (E 21 ) = D\ = D x , (126) 
Hi = 2N X - n , (12c) 
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where we have denoted z \ = x } Nl = N x . Note that we have obtained the well known 
[21] realization of the U q (sl(2)) representations with X + = T 2 (Ei 2 ), X~ = ^(£21), 
H = Hf, depending on the representation parameter r\, (ro being cancelled as expected). 
For q = 1 this coincides with (2) setting H = 2 Jo, X^ 1 = J±. 

Next we take n = 3 setting Nf = Ni, D 3 = Di, z\ = z, z\ = y, r = r 2 = r\ + r 2 ■ 
Besides this we renormalize the generators r 3 (Ei 3 ) and Y^(E^i) so that they obey the 
standard relations [22], [20] (these are different in [15], cf. (6.17), (6.20)) : 

r 3 (£ , i 3 ) = r 3 (£ , i2)r 3 (£'23) — q 1/2 t 3 (e 23 )t 3 (e 12 ) , (i3a) 
t 3 (e 3 i) = t 3 (e 32 )t 3 (e 21 ) — q~ 1 / 2 r 3 (£ 2 i)r 3 (£ 32 ) . (136) 

Thus, we have: 

Ts(E 12 ) = V 2 (E 12 ) q k( N 't-Nt) + ^^fe-Bn) Z 3 D * = 

= x[n -N x ] qi( N *- N v) + zD y q^- 2N ^ , (14a) 
T 3 (E 21 ) = T 2 (E 21 ) q^-Ni) + q ^ 2 (E 22 -E ll} z s D s = 

= D x q *( N *- N v) + yD z q^~ 2N ^ , (146) 
Hi = 2N 2 - n + Nl - Nl = 2N X - n + N z - N y , (14c) 
Hi = 2Nl -r 2 + Nl - Nl = 2N y - r 2 + N z - N x , (14a 7 ) 

r 3 (£ 31 ) = q WT 3 (E 22 )-T 2 (E 22 )) D 3 = Dzq l 4 ( ri + r -N x + 2N y ) ^ ^ } 

T 3 (E 32 ) = qF'(En) B\ = D y q^ ro + N ^ , (14e") 



t 3 (e 13 ) = q^(^)-w 22)) ,3[r 2 ( jB33 )-r2(£ 11 )-^iv, 3 ] - 

k=l 

- q-^ T3(E ^q ai2 4^2(E 12 ) = 

= z[r - N x - N z - N y ]q^ N *- r '- 2N y- ro) - 

- xy[ ri - Nx ]q^ 2r i- r ° + N *- N z- 3N y +1 "> (14/') 

2 

r 3 (£ 23 ) = q-^ 2(E ^ z 3 2 [r 2 (E 33 )-r 2 (E 22 )-Y,N 3 k ] - q a21 zlv 2 (E 21 ) = 

k=l 

= y[r 2 +N x -N z - N y ] q-^+ro) _ 

_ z £f x q -$(2r- ri + r + N x -N z -Ny+l) _ (14/") 

We now rescale the generators E 3 i 7 Ei 3 (i = 1,2) so as to absorb the parameter ro- (Such 
a rescaling should be done also for the general U q {gl{n)) case.) Thus the realization of 
U q (sl(3)) depends only on the parameters ri 7 r 2} as in the classical case. The latter is 
obtained from (14) by setting q = 1, (d a = '■ 

T 3 (E 12 ) = x( ri -xd x ) + zd y , (15a) 



y 

r 3 (£ 2 i) = d x + yd z , (156) 

Hi = 2xd x - yd y + zd z - r 1 , (15c) 

Hi = -xd x + 2yd y + zd z - r 2 , (15d) 

r 3 (£ 3 i) = d z , (15e') 

r 3 (£ 32 ) = d y , (15e") 

T 3 (E 13 ) = z(r - xd x - zd z - ydy) - yx(r! - xd x ) (15/') 

r 3 (E 23 ) = y(r 2 + xd x - zd z - ydy) - zd x . (15/") 



Of course, (15) can be derived from the induced representation picture [19]. 



3. Reducibility of the representations and invariant subspaces. 

3.1. Let us apply the realization (6) to the function 1. Using the fact that 

N? 1 = = £>f 1 we have: 

i-l 

T n (Eu)l = Y, r J = r * _1 > H i 1 = " r * ' l ^ n > ( 16a ) 

j = 

Y n (E ni ) 1 = , % < n , (166) 

T n (Eij)l = Tn-^Eij) 1 = ••• = T,^) 1 = 0, j<i<n, (16c) 

r„(E ij ) 1 = Tn-^Eij) ! = ••• = r,^) 1 , i < j < n , (lQd) 



where 



r n (^ m ) i =^i^ =i+1 r 2^ =1 r ) z n \ (r . i + ... + 



r n -i 



n-1 

V a? ' < ^s(E ia ) 1, i<n 

s=i+l 



= gHE::: +1 '- 4 - i -E;: 1 1 '- 4 - 1 ) ^('■ b - i -e::>* b )+i 



(16e) 



(17) 



It is straightforward to obtain the explicit expressions for T n (Eij) 1, i < j < n, applying 
recursively (16c/, e). In particular, we have: 

r„(£ iii+ i) 1 = T i+1 (E i>i+1 ) 1 = 9 -^El=i r * _1 [ ri ] 2r «+ 1 . (18) 



Thus, we have obtained a lowest weight module (LWM) with lowest weight vector 1 , 
(it is annihilated by the lowering generators T n (Eij) } j < i < n), and lowest weight A 
such that A(Hi) = —ri } (cf. (16a)). Generically this LWM is irreducible and then it is 
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isomorphic to the Verma module with this lowest weight. The states in it correspond to 
the monomials of the Poincare-Birkhoff-Witt basis of U q (Q + ), where Q + is the subalgebra 
of the raising generators. This is isomorphic to the monomials in the variables z k . When 
the representation parameters V{ or certain combinations thereof are nonnegative integers 
our representations are reducible. The full analysis of this is postponed to the sequel of 
this paper. Below we consider in detail the cases n = 2 and n = 3. 

3.2. We start with n = 2 (though this example is well known). Let us apply 

H = T 2 (E 11 - E 22 ), X+ = T 2 (E 12 ), X~ = T 2 (E 21 ) to the function 1. We use 
the fact that N x 1 = = D x 1. Thus: 

H 1 = -r , X+ 1 = x[r] , X~ 1 = . (19) 

Thus, we obtain a lowest weight module with lowest weight vector 1 and lowest weight 
A such that A(H) = —r. All states are given by powers of x , i.e., the basis is x k with 
k G -2v+ and the representation is infinite dimensional. The action of U q {s!(2)) is given by: 

X+x k = [r-k]x k+1 , X~x k = [k^- 1 , Hx k = (2k-r)x k . (20) 

Clearly, if r ^ this representation is irreducible. Furthermore all states may be 

obtained by the application of X + to the LWV, i.e., : 

(X + ) k 1 = x k [r][r -l]...[r - k + 1] , k G 25+ . (21) 

Let r G -2v+ , then (X+) r+1 1 = X+x r [r]\ = 0. Thus, the states x k with 
k = 0, 1, ... r form a finite-dimensional subrepresentation with dim = r + 1. Note that 
the complement of this subrepresentation, i.e., the states x k with k > r, is not an invariant 
subspace. 

Clearly, any polynomial in if, X^ 1 , will preserve this invariant subspace and thus would 
be a quasi-exactly solvable operator. 

The invariant subspace may be obtained as the solution of either one the following 
equations: 

(X+Y+ 1 f(x) = , (22a) 
(X-y +1 f(x) = , (226) 

in the space of formal power series f(x) = fJ>kX k ■ Note, however, that only (226) 

(which is enough) was expected - this is an artefact of n = 2 simplifications. Indeed, only 
the operator in (226) has the intertwining property (as in the classical case [19]): 



(X-) r+1 T 2 (X) r = T 2 (X) r , (X-) r+1 , r ' = -r-2, 



(23) 



where X = X 1 * 1 , and T 2 (X) r is from (12) with explicit notation for the representation 
parameter of the two representations which are intertwined. 

3.3. To the end of Section 3. we discuss U q (sl(3)). 
3.3.1. Let us apply (14) to the function 1 : 



r 3 (£i 2 ) i 


r i 

= x[ri\ , 


(24a) 


r 3 (£ 2 i) i 


= o , 


(246) 


#i i 


= ~ri , 


(24c) 


H 2 1 


= -r-2 , 


(24d) 


r 3 (£ 3 i) i 


= o , 


(24e') 


r 3 (-E 32 ) i 


= o , 


(24e") 


r 3 (£i 3 ) i 


= q~^ ri z[r] — q^ ( - 2r2+1 "'yx[ri] 


(24/') 


r 3 (-E 23 ) i 


= y[r-2\ , 


(24/") 



Thus, we obtain a lowest weight module with lowest weight vector 1 and lowest weight 
A such that A(Hk) = — All states are given by powers of x } y } z , i.e., the basis is 
generated by x J z k y £ with j, k,£ £ The action of U q (sl(3 j) is given by: 

T 3 (E 12 ) xH k y £ = [n-j] q^ k ~ £) x 1+1 z k y £ + [£]q^~ 2 ^ xH k ^y £ ~ X , (25a) 



r 3 (£ 21 ) x j z k y e = [j] q^ k ~ £) x^ 1 z k y £ + [k}q^~ 2 ^ x^z k - 1 y £+1 , (256) 

H! x J z k y £ = (-n + 2j - t + k) x J z k y £ , (25c) 

H 2 x J z k y £ = (-r 2 - j + 2£ + k) x J z k y £ , (25d) 

T 3 (E 31 ) xH k y £ = [k]q^-^ 2£ ^ x^z k - 1 y £ , (25e') 

r 3 (E 32 ) x^z k y £ = [i]q& x^z k y £ - 1 , (25e") 

r 3 (E 13 ) x j z k y £ = q^- r '- 2£) [r - j - k - £] x J z k+1 y £ - 

_ 9 *(2r 2 +j-*-3*+i) [ri _ ;] x Hi z k y t+i ^ (25//) 

r 3 (E 23 ) x j z k y e = q~^ [r 2 + j - k - £] x J z k y £+1 - 



In this Section we show the following results which parallel the classical situation (cf. [19]): 

1. If ri, or r 2 , or r+1 £ .2v_|_ this representation is reducible. It contains an irreducible 
subrepresentation which is infinite-dimensional, except when both r\ , r 2 £ ; 

2. If ri , r 2 , r + 1 ^ this representation is irreducible and infinite dimensional. 
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3.3.2. Clearly, if r\ £ the representation (25) becomes reducible : the monomials 
x^z k y e with j < ri form an invariant subspace since from (25a,/') we have: 

r 3 (£i 2 ) x ri z k y e = [£]q~^ x^z^y'- 1 , (26) 

r 3 (£i 3 ) x ri z k y e = [r 2 -k-£]q-^ e x r ^z k+1 y e , (27) 

and all other operators are either lowering or preserving the powers of x. This invariant 
subspace may be described as the solution of the following q - difference equation: 

(D x r +1 f(x,y,z) = 0. (28) 

Note that the operator in (28) has the intertwining property (as in the classical case [19]): 

(D x ) r * +1 T 3 (X) ri , r2 = r 3 (X) r , ;r , (D x ) r * +1 , r[ = -n - 2 , r' 2 = r + 1 , (29) 

where X = En — Ei+i^+i, Eij } i ^ j, r 3 (X) ri>r2 is taken from (14) with explicit depend- 
ence of the representation parameters of the two representations which are intertwined. 

The subrepresentation obtained is infinite dimensional if r 2 -2v+ since the powers 
of y,z are still unrestricted by (25/',/"). 

3.3.3. If r 2 £ -2v+ the representation in (25) becomes reducible. In the classical case 
(q = 1) the equation which singles out the invariant subspace is [19]: 

(xd z + d y ) r2 + 1 f{x,y,z) = 0, q = l. (30) 

For the quantum case we have the following expression : 

g D 2 r2+1 f(x,y,z) = 0, (31a) 

k , 

q v? +1 = ]T x k ~ s D k -° d; q k^^-^)H(s-k)N y -kkN x ^ (316) 

which coincides with (30) for q = 1. The invariant subspace is infinite dimensional if 

n i 2Z+. 

As in the classical case [19] the explicit form (316) of this operator may be checked 
by the intertwining property: 

q V? +1 T 3 (X) ri , r2 = r 3 (X) r[ y 2 q V?+\ r[=r + \, r 2 = -r 2 -2, (32) 

where X = E tJ} i ^ j } E tt - E i+lji+1 , r 3 (X) ri>r2 is from (14) as in (29). 
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3.4. In this subsection we consider the case when both r& £ . 

3.4.1. Let rk £ -2v+ , k = 1,2. Then there is a finite dimensional irreducible subspace 
of dimension: 

d ri ,r 2 = ^(ri + l)(r 2 + l)(r + 2) . (33) 

Thus, we recover the complete list of the finite dimensional irreps of U q (sl(3)) and 
SX(3), and by default, also the complete list of the finite dimensional unitary irreps of 
U q (su(3)) and SU(3) (we have assumed that q is not a nontrivial root of 1). Below we 
give explicitly the basis of this subspace. 

3.4.2. Let us consider first the special case r\ = 0, r 2 = r £ Z+. (This is the case given 
in [17] for q = 1.) Then there is no x dependence in the irreducible subrepresentation as 
we argued above. Since ^(£12) 1 = 0, all states can be built in the following way. Apply 
first r3(£^23) to the function 1: 

(T 3 (E 23 )Y 1 = [r][r-l]...[r-^ + l] y e (34) 

thus we obtain r + 1 states (^3(^23 )) r+1 1 = 0). Then we apply the operator ^(£13) 
to each of these states: 

(r 3 (£i3))V = q~? ki [r-£][r-£-l]...[r-£-k + l] y e z k (35) 

thus for each fixed £ we obtain r + 1 — £ states. 
Thus the number of states is (cf. (33)) : 

r r— I ^ 

EE 1 = 2 {r + 1)(r + 2) = do ' r - (36) 

£=0 k=0 

3.4.3. In a similar way we obtain the case r\ > 0, however, hypergeometric finctions 
enter the game. Apply first ^(£12) to the function 1: 

(r 3 (£i 2 )) J 1 = [r 1 ][r 1 -l]...[r 1 -j + l] X i (37) 

thus we obtain r\ + 1 states ((r3(£i2)) ri+1 1 = 0). Then we apply the operator ^(£13) 
to each of these states: 

k 

(r 3 (£i 3 ))V = i- 1 ) 3 ^ V • (38) 

s=0 

Applying ^(£13) once more we obtain a recursion relation: 
Pk+i,s = q^- ri - s) [r-j-s-k]/3 ks + q ^+j-k-s+2) ^ _ ■ _ $ + (39) 
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with the convention that @k,k+i = ftk,-i = 0. Solving this, we have : 

[ri - j]l[r - J - s]l 



P 



k,s 



qkHj-rx) + \s{2r- ri +2-k) 



[ri - j - s]\[r - j - k]\ 



Note that: 



, if 



> min(A;, r\ — j) for r ^ r\ 



> k 



for r = r\ 



Thus we have: 



(T 3 (E 13 )) k T 3 (E 12 )) j 1 



x J z k q \ k (j-ri) 



[nV- 



[r-J-k}\ ^ 



l r ~ ■? ~ 3 J- i S (2r-n+2-*) 

[ri — j — s]l 



z 



(40) 



(41) 



x J z q< 



[r-j- &]![n - ]}\ 
x 2F«(-k,j 



ri ;j-r; q ^-r 1+ 2-k) ^ 



for r ^ r\ 



(42) 



x 3 z kq\k(j r) 



-j-k]l 



1 



9 



±(r+2-fc) £1 



for 



r = n 



where 2-^i 9 is the g-hypergeometric function, (1 + w) k is the g-binomial. From (42) we 



have: 



(r 3 (£i 3 ))V =0, k>r-j . 



(43) 



Thus, we have at this stage (ri + l)(2r + 2 — ri)/2 states. Now it remains to apply 
(r 3 (E 2 3)) £ to (r 3 (£i 3 )) fc r 3 (£i 2 )) J 1 to obtain the remaining r 2 (r 1 + l)(r + l)/2 states. 

Let us denote v^j = T 3 (Ei 3 j) k r 3 (.Ei2)) J 1. If r = r\, (r 2 = 0), we have: 
T 3 ( J B 23 ) v k j = 



[j> J " 1 z fc+1 g^ ((fc " 1)(r+1) " j(fc+1)) 



j-k]\ 



1 — q 4 



(r+l-k) X V 



k+1 



(44) 



i.e., r 3 (£^23) just transforms the states obtained so far. Thus, the number of the states is : 



EE 1 = i(r + l)(r + 2) = d r ,o 



(45) 
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as expected (cf. (33)). 

Statement. This subspace may be characterized as the solution space of two equations, 
namely, (28) and: 

(xD z q-^ N y + N * ) + D y q^ N *- N *- r) ^ f{x,y,z) = 0, (46) 
which is (31) for r 2 = 0. 

Proof. Let us write f(x,y,z) = Y2j k ie% + ^jki^^ z k y l and apply (46). Below we 
consider Xjki for j, /.-.'• '//. with the convention that Xjki = if j < 0, or k < 0, or £ < 0. 
We obtain the following recursion relation for Xjki- 

[k + l]X^ hk+ i,i + [e + l]\ jk ,i+iq* (k+e - r) = 0. (47) 

We note that from (28) follows that Xjki = if j > r\ = r. From this and (47) follows 
that Xjki = if j ' + k > r, (since Xjki ~ ^r+i,A;+j-r-i,^-j+r+i = with proportionality 
coefficient being non-zero and finite for the hypothesis in consideration). Next we note 
that Xjki = if £ > j, (since Xjki ~ X-i^+j+i^-j-i = 0). We also note that Xjki 
are proportinal to each other for k + £ fixed. In order to ensure the above properties 
automatically we are prompted to make the substitution: 

j — > j -\- s } k — >■ k — s , £ — > s } with s £ , s < k . (48) 

Solving (47) with this we get: 

X, ■ .-./, ,,, = (-iy q^^X^o , (49) 

where Xjko are only subject to the vanishing condition for j -\-k > r. The latter is achieved 
if we suppose that Xjko ~ l/[ r— j~ k] q \ ■ Thus we see that up to normalization A j+ Sj k-s, s 
coincide with fiks for r = r\. • 

3.4.4. Thus, it remains to treat the case when r& £ IV. First we note another case 
where the r3(£^23) transforms the states: 

T 3 (E 23 ) Vr-jj = -q~^ r+1) [j] Vr-j+u-! , (50) 



ID 



Next we shall use the following general formula valid for arbitrary r& 
v ekj = T 3 (E 23 f T 3 (E 13 ) k T 3 (E 12 ) j 1 = 



ES->r (!) (I) ,* 



( 4{(j-ri)A:-^+(s-n)(ri + 2r 2 -A:-^+2)} 



r g (n + l)r g (r-j-s + l) r g (r 2 +j +s - At -n + 1) [j + 



r g (ri -j -s + l)r g (r-j -Ar + 1) T q (r 2 + j + s - k - £ + 1) [j + s - n}\ 

X X J+S_n y i+s ~ n 

(51) 

Note that for r 2 G -2v the ratio r g (r 2 + j + ,s — At — n + l)/r g (r2 + j + .s — At — ^ + 1) means 
just [r 2 + j + s — At — n] [r 2 + j + s — At — n — 1] . . . [r 2 + j + s — At — / + 1] . 

A basis for our representation space is given by vgkj iff ^ + At + j ' < r, < j < r\ , 
< £ < r 2 . To show this it is enough to demonstrate that vgkj with any of the conditions 
violated is a linear combination of the representation states (which combination may be 
zero) - cf. e.g., (44) and (50). In particular, by induction, it is enough to show this for 
the boundary cases: £ = r 2 -\-l } £-\-j-\-k = r-\-l. We give explicitly the boundary case 
£ = r 2 + 1. Let us introduce the notation 

A = T 3 (E 23 ), B = T 3 (E 12 ), C = T 3 (E 13 ). (52) 

Using the commutation relations it is easy to show that: 

A'C f ' = , C h B j = q^B J C k , (53) 

BiAt = 1 t(t) ( J t) m**"^*-*^-* , A = min{^}. (54) 
t=o ^ 7 ? ^ 7 ? 

The proof of the result stated above now follows very easily from (53), (54) and the fact 
that A r2+1 1 = 0. In fact 

7 ' 02+ 1 )* 7 

A r 2 + l C k B J x = q —^~ C k X 

x ( Si ^ +1 -E( r % +1 ) (i) m^^^-'c'B^ 1, 

(55) 

here A = min{r 2 + 1, j}- The first term in the parenthesis gives zero when applied to 1, 
while the second term yields a linear combination of states of the required form (£ < r 2 ). 
This ends the proof. The boundary case £-\-j-\-k = r-\-l may be treated in a similar way. 
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3.5. Next, we discuss the case when r + 1 £ , but r& ^ .2v_|_. In the classical 
case there is a nontrivial invariant subspace singled out by the following equation [19]: 

V r + 2 f(x,y,z) = 0, (56a) 

t>; +2 = e , ( ;] r _ , ( r + 2 ) (xtr 2 ~ s (xtr 2 A + r , (566) 

X+ = <9 X , X+ = xd z + a, . (56c) 

For this equation one uses the singular vector in (8.41) of [19]. The necessary singular 
vector in the q case is obtained from the classical case by just replacing the numbers with 
g-numbers [20]. However, in order to generalize this to the q - case we first rewrite (566) 
with the explicit substitution of (56c). We obtain: 

v ' +2 =-% r(-,f+'+ 1 2'-l ) ) (T) d ' W'+W-' 9 '^' = ( 57o » 

(r + 2)! 2 T(-1 -n) 



EE 



s\t\(r + 2 - s - t)\Y(r - n + 2 - s) 

r+2-s-t 

d r + 2 - f d\ d f x Yl ( xd * -< + !-«)• (576) 



u = l 



Now we have the following expression for the q - difference equation in our case: 

q V r + 2 f(x,y,z) = 0, (58a) 



r+2 v^v^ g^ r+2 -*- 2s ) ri+ iK r+2 )* + iK r+1 ^- 4 ) r g (-i - ri) 

— X 

(586) 



,v 3 

s =o t=o 



9 3 Z^Z^ [s]![<]![r + 2-s-<]! r g (r-ri + 2 - s) 



r+2-s-t 



D r z +2 -< D\D\ J] [N x - t + 1 - U] q ^^N x + LN v + ^2l Nz 



u = l 



As in the classical case [19] the explicit form of this operator may be checked by the 
intertwining property: 



q v r + 2 r 3 (x) ri , r2 = r 3 (x) ri , r / q v r + 2 , r [ = -r 2 - 2, r ' 2 = -n - 2 , (59) 

where X = E tt - E i+lji+1 , E tJ} i ^ j, r 3 (X) ri;r2 is from (14) as in (29), (32). 

The states in the subrepresentation are given by vgkj, with £ } k } j £ -Zv+, k < r + 1. 

Let us illustrate the above by the simplest classical example ofr = — l,g = l. The 
equation (56a) has the form: 

((ri -xd x ) d z - d y d x ^j f(x,y,z) = . (60) 
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It is easy to see that all solutions of (60) are given by: 

f jt (x,y,z) =x>y e g n! (j - n)\ (I - n)\ T(ri - j + 1 + n) (^) = (fJ1) 

= x> y l 2 F x {- h -i- n - j + 1 ; — ) , j,£eZ+ . 

xy 

The same states are valid in the g-case since one has from (51): 



Vl0j\r=-1 = 2^ \ n ) q 



r,(ri + i) r g (i + n-j+i) xJ _ n zn y£ _ n 



T q (n -j + 1) r g (l + n -j + n) [j -n]\ (62) 

T q ( ri -j + l) r q (j- ri -e) 



x j y l 2 F«(-j,-£ ; n - j + 1 ; ^ (ri+ ^ — ) . 

xy 

Alternatively one may check that this is the general solution of (58) for r = — 1. 



3.6. Finally, we discuss the case when ri 7 r 2 G -2v , r + 1 = ri + r 2 + 1 G -2v+ , 
rir 2 G — -2v+. 

We first take the case r\ G -2v+ , thus r 2 G —IV and r < r\. Then we restrict r 2 < — 1. 
In this situation there is an invariant subspace which is singled out by the two equations 
(28) and (58), however, in (586) the two factorials T q ( — 1 — ri)/Y q (r — r\ + 2 — s) are 
replaced by the regular expression: ( — l) r+1 ~ s [ri — r — 1 + ■s]!/[ r i + 1]' • The states 
in the subrepresentation are given by vgkj, with £ } k } j G -2v+, and either k } £ > 0, 
—r 2 — 1 < j < ri, or ^ > 0, < At < r + 1, < j < — r 2 — 2. (The last statement 
is a conjecture proved in partial cases.) 

Then we take the case r\ G -2v+ , r 2 = — 1. In this situation already the classical 
operator (56) has to be modified. It has the form (up to overall multiplicative non-zero 
constant) [19]: 

V r + 2 = (X+Y+ 2 (X+Y+ 2 , (63) 

which can be obtained from (56) by multiplying (566) with r 2 + 1 and then taking the 
limit r 2 — > —1. Thus, it is a composition of the two operators (28) and (31). Since these 
operators are invariant, the invariant subspace is given by the solutions of (28) only. The 
latter conclusion may be reached also from the previous subcase taking r 2 — > —1; then the 
second set becomes empty and we are left with what is claimed here. 
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We then take the case r 2 G -2v+ , v\ < — 1. In this situation there is an invariant 
subspace which is singled out by the two equations (31) and (58), The states in the 
subrepresentation are given by vgkj, with £ } k } j G -2v+, and either k } j > 0, —r\ — 1 < 
I < r 2} or j > 0, 0<k<r-\-l } 0<£< —r\ — 2. (The last statement is a conjecture 
proved in partial cases.) 

Finally we take the case r 2 G -2v+ , r\ = — 1. In this situation again the classical 
operator (56) has to be modified. It has the form [19]: 

v;+ 2 = (x+y+ 2 (x+y+ 2 , (64) 

which can be obtained from (56) by multiplying (566) with r\ + 1 and then taking the 
limit r\ — > —1. Thus, it is a composition of the two operators (28) and (31). Since these 
operators are invariant, the invariant subspace is given by the solutions of (31) only. The 
latter conclusion may be reached also from the previous subcase taking r\ — > —1. 

Clearly, if ri,r2,r + 1 (ji the representation in (25) is infinite dimensional and 

irreducible since by (25a) the powers of x are arbitrary, by (25/") the powers of y are 
arbitrary, by (25/') the powers of z } x } y are arbitrary. 



4. Newton diagrams. 

It this Section we give a visualization of the representation spaces. Each state is 
represented by a point on an integer lattice in n(n — l)/2 dimensions, i.e., on /Z^ n 1 ^ 2 . 
For a finite-dimensional subrepresentation the number of these points is finite and the hull 
of these points is a convex polyhedron in IE!^ n 1 ^ 2 . Such a polyhedron (not necessarily 
convex) is called a Newton diagram [18]. In the present context this notion was introduced 
in [17], where also some examples in the case of functions in one and two variables were 
given (for q = 1), when the figures are planar (polygons). Below, we give explicitly 
the Newton diagrams for n = 3. Moreover, we introduce also infinite Newton diagrams to 
depict the infinite-dimensional nontrivial subrepresentations. 

4.1. Finite Newton diagrams for n=3. Fix r^ G -2v+ . Then the Newton diagram 
is given by the points with integer coordinates j,£,k in such that: 

< j + k + £ < r, (65a) 

< j < ri , (656) 

< t < r 2 , (65c) 

cf. the Figure. The polyhedron formed by these points is planar only for r\ = or r 2 = 
in which case it is a triangle (only (65a) is relevant since r = r 2 or r = r\). (The case 
r\ = was given in [17].) 
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Fix a point j,£, k. This is represented by the state vgkj- Then, the number of states 

is: 

n r—j min(r— k— j,r2 ) r\ r\—j r 2 r\ r—j r—k—j 

EE E 1 = E E E 1 + E E E 1 = 

j=0 k=0 £=0 j=0 k=0 £=0 j=0 k=r 1 -j+l £=0 



(66) 



(ri + l)(n +2)(r 2 + 1) (n + l)r 2 (r 2 + 1) _ 



as expected (cf. (33)). 

Note that such diagrams have an advantage over the usual weight diagrams for sl(3) 
and su(3) which are degenerate. For instance, consider the adjoint representation obtained 
for n = r2 = 1. The weight diagram consists of two orbits of the Weyl group, one with 
six points with multiplicity one, and the other with one point with multiplicity two. To 
the latter point in our diagram correspond the two states: 

v wl = ([2] q xy - q~ x z) , (67a) 

^oio = ([2] q z - qxy) , (676) 



which are linearly independent. 



4.2. Infinite Newton diagrams for 11= 3. Here either r\ or r 2 ^ -2v+ and the 

considerations run in parallel with Subsections 3.2, 3.4, 3.5. Below j } £ } k £ 

4.2.1. For r\ £ and r 2 , r + 1 ^ the Newton diagram is given by the points with 
coordinates: 



< k , 

< j < r 1 

< i . 



(68a) 
(686) 
(68c) 



4.2.2. For r 2 £ and ri, r + 1 ^ ^_|_ the Newton diagram is given by the points with 
coordinates: 

< k , (69a) 
< j , (696) 
< ^ < r 2 . (69c) 
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4.2.3. For r + 1 £ .2v_|_ and ri, r 2 ^ -2v+ the Newton diagram is given by the points with 
coordinates: 

< k < r + 1 , (70a) 

< j , (706) 

< ^ . (70c) 

4.2.4. For ri,r + 1 £ ^+ and r 2 + 1 £ —IV the Newton diagram is given by two sets 
of points with coordinates: 

< k , (71a) 

- r 2 - 1 < j < n , (716) 
< t . (71c) 
< k < r + 1 , (72a) 
< j < -r 2 - 2 , (726) 
< t . (72c) 

4.2.5. For r\ = r + 1 £ ^_|_ and r 2 = —1 the Newton diagram is given by (68). It 
can be obtained formally from the previous case by setting r 2 = —1 , then (71) coincides 
with (68), while (72) is empty. 

4.2.6. For r 2 ,r + 1 £ and r\ + 1 £ —IV the Newton diagram is given by two sets 
of points with coordinates: 

< k , (73a) 

< j (736) 

- r 1 - 1 < j < r 2 , . (73c) 
< k < r + 1 , (74a) 
< t , (746) 
< j < -r 2 - 2 . (74c) 



4.2.7. For r 2 = r + 1 £ ^_|_ and ri = —1 the Newton diagram is given by (69). It 
can be obtained formally from the previous case by setting n = — 1 , then (73) coincides 
with (69), while (74) is empty. 
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Figure Caption. 

Newton diagram for the finite-dimensional representations of U g (sl(3)). 
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Newton diagram for the finite-dimensional representations of U g (sl(3)). 



